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The rate of metastable decay in nonequilibrium systems is expected to display scaling behavior: 
i.e., the logarithm of the decay rate should scale as a power of the distance to a bifurcation point 
where the metastable state disappears. Recently such behavior was observed and some of the earlier 
predicted exponents were found in experiments on several types of systems described by a model of 
a modulated oscillator. Here we establish the range where different scaling behavior is displayed and 
show how the crossover between different types of scaling occurs. The analysis is done for a nonlinear 
oscillator with two coexisting stable states of forced vibrations. Our numerical calculations, based 
on the the instanton method allow the mapping of the entire parameter range of bi-stability. We 
find the regions where the scaling exponents are 1 or 3/2, depending on the damping. The exponent 
3/2 is found to extend much further from the bifurcation then were it would be expected to hold as 
a result of an over-damped soft mode. We also uncover a new scaling behavior with exponent of ~ 
1.3 which extends, numerically, beyond the close vicinity of the bifurcation point. 

PACS numbers: 



I. INTRODUCTION 

The problem of noise-induced escape from a 
metastable state has been studied in various contexts, 
from nucleation to chemical reactions to switching in 
nanomagnets. The escape rate often has the Arrhenius 
form W oc exp(—R/D), where D is the noise intensity 
(temperature, in the case of thermal noise). Of central 
interest both for theory and experiment is the activation 
energy R. Starting with the Kramers paper [l[ much 
work has been put into calculating R and the prefactor 
in the escape rate for various physical systems. From 
this point of view, it is particularly important to reveal 
generic features of the activation energy R, such as scal- 
ing behavior with the system parameters. The onset of 
the scaling behavior was noticed first for systems close 
to thermal equilibrium @, 0] . Scaling occurs also for sys- 
tems far from thermal equilibrium, which generally do 
not have detailed balance 

[1 i, 0, S H • Nonequilibrium 
systems may display several types of scaling behaviors 
and crossovers between different scalings with varying 
system parameters. 

A nonequilibrium system that has attracted much at- 
tention recently is a Duffing oscillator modulated by a 
periodic force 

q + 2Tq + u\q + ~/q 3 = A cos (u F t) + f(t) (1) 

with high quality factor Q = |p- The term / is white 
noise with property (/(*)/(*')) = 2BS(t - t'). As a re- 
sult of the modulation the oscillator may have two or 
more coexisting states of forced vibrations. Noise leads 
to switching between these states. In this paper we will 
use the following two parameters: dimensionless friction, 



CI 1 = T /\u>f — u)q\ and a dimensionless driving strength, 

VP = Ay/3j/32w F \w F - coo 
f2 _1 , the region of bi-stability is closed (Fig. [TJ . 
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FIG. 1: Region of bi-stability in the parameter space of the 
driving, */j3 and friction, Cl^ 1 . 

Upon varying \/p at a fixed CI one would encounter an 
amplitude response curve such as that depicted in Fig. [2j 
which exhibits a region of bi-stability between two bifur- 
cation points \J (3 l B and \J~P%- It is critically important 

to note that at any finite CI , bi-stability will exist only 
at a non-zero (3. Therefore the problem of transition 
rates between two stable attractors is a non- equilibrium 
problem! 

Several types of scaling behavior of the switching acti- 
vation energy have been seen for such different systems 
as modulated nanoelectromechanical (NEMS) [lfj, 11, 
[l6j . and microelectromechanical (MEMS) oscillators 13 1 
(various sources of noise in the nanomechanical systems 
have been thoroughly described in M), Josephson junc- 
tions [3], optically trapped atoms [la ], and supercon- 
ducting resonators [TtJ - all of these systems are modelled 
as Duffing oscillators. Not only is this scaling interesting 
on its own, but it also underlies various applications, an 
important example being quantum measurements [14j . 
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FIG. 2: Amplitude versus \ffi for some given fT 1 / 0. The 
upper and lower branches corresponds to a stable steady state, 
while the middle branch corresponds to a steady state which 
is unstable with respect to small perturbations. Similar curve 
exists for phase versus y 7 /?. A bifurcation value of ^fp will 
be denoted by \/(3b; there are two of these: \J /3g and \f f3 l B . 
Noise-induced transitions from low amplitude to high ampli- 
tude branch will be called "down — > up transitions" , and those 
from high amplitude to low amplitude branch, will be called 
" up — > down transitions" . 



Let f] = \\/~j3 — V/Sb] (the letter 77 is reserved for de- 
noting (3 — Pb in recent work by one of the authors, see 
[26| for example). In certain cases, the switching acti- 
vation energy obeys a power law: R ~ fj^. There exist 
regimes when the escape problem can be reduced to that 
of escape over a 1-dimensional potential barrier. Analyt- 
ical approximations are valid in two such regimes - close 
to bifurcation points, where £ = 3/2 and at at low Q , 
where £ = 1 are summarized in Section II. A. A new geo- 
metric explanation of these £ = 3/2 and £ = 1 exponents 
has recently been given (2?]]. These authors make use of 
a remarkable feature of this system - that the bifurcation 

' f3 l B is " non-local" while the bifurcation at J P B is 



"local", as will be explained in Section II. B. 

In this paper we present a detailed theoretical study of 
the scaling and scaling crossovers of the activation energy 
for a Duffing oscillator. Using numerical techniques, we 
map out the entire bi-stability region shown in Fig. |TJ We 
identify part of this bi-stability region where the £ = 3/2 
and £ = 1 scaling laws hold, where there are crossovers 
between these scaling laws, as well as regions of other 
scaling behaviors not previously mentioned in the litera- 
ture. 



II. BRIEF REVIEW OF EXISTING THEORY 

We will analyze the system in the canon- 

ically transformed phase space (Q,P), with 

q = 2 3 ^ 2 loo ^J^f^(Q cos (u F t)+P sin {u F t)) and 



= 2 3 / 2 w F Wo 



3 7 



(— Q sin (cup t) + P cos (copt)). In 



q 

the underdamped limit when Q » 1, the motion of q is 
approximately simple-harmonic, with slow "envelope" 
evolution represented by P and Q. The resulting 



equations for P and Q are approximately independent 
of time - the fast time-dependent corrections average 
to zero to lowest order in Q^ 1 . In terms of slow time 
T = QT x t (subsequently, the dot denoting d/dT), it 
follows the following evolution [28| : 



Q = K Q (Q,P)+n Q (T) 
P = Kp(Q,P) + n P (T) 



(2) 
(3) 



where K Q {Q,P) = J£ 



Cl^Q and K P {Q,P) 



.SR. 

dQ 



ft 1 P, g(Q, P) is a Hamiltonian given by g = \{Q 2 +P 2 — 
l) 2 — ^/]3Q and fi" 1 is an effective damping which tends to 
arrest the amplitude to zero. The n^s are stochastic terms 
which will be modelled as white noise: (rii(T)nj(T')) = 

2D5 lJ S(T - T) where D = ( JjF F 2 ) ^~ 2 - The terms 

inside (...) have been defined on the first page and are all 
parameters of the oscillator. In the absence of noise, an 
initial condition (Qo, Po) will evolve into one of the two 
attracting fixed points (FP), or " attractors" . The third 
FP is of saddle type - its amplitude forms the middle 
branch in Fig. O The separatrix between the basins of 
attraction of each of the two attractors is a set of initial 
conditions that flows into the saddle FP via its attractive 
eigen-direction. 

In the presence of noise, an initial condition 
(Qo,Pq) follows a continuous random walk in the 
deterministic flow K. The probability density 
(PDF) p(x, y) of reaching (Q, P) from some ini- 
tial (Qo,Po) satisfies the 
(FPE): §§ - {g, P } = V • 
solution can be written as a functional integral of 



Fokker-Plank equation 



rp - 



DVp) . The 



exp 



iJS; p p\) c (Q'( T )> p '( T )) dT ° ver a11 p aths 



D J(Qo,Po) 

H, where C = 



P-Ki 



One 



can re-express all paths as variations around the optimal 
path that minimizes the exponent. Then p(Q, P) = 

p'(Q, P) exp [-iming, C (Q'(t), P'{T)) dT = 

p'(Q,P)exp -S(Q,P)/D] where the prefactor p'(Q,P) 
comes from performing path integrals over variations 
around the optimal path. This pre-factor is nonexponen- 
tial because the path integral over variations is to the 
lowest order a Gaussian path integral. We are interested 
in the exponential factor, because in the limit of weak 
noise it will dominate transition rates. 

Now, the probability of undergoing a large fluctuation 
from a FP to cross the separatrix is proportional to the 
probability of being found in the most likely place where 
the exit traffic takes place - the saddle point (Q S ,P S ) 
flU - Therefore, the rate of escape W is to be approxi- 
mated by p(Q s ,P s ). To calculate the probability of es- 
cape up to this accuracy, we need in principle to solve 
Euler-Lagrange equations to find a path that leads from 
the attractor to the saddle and evaluate the action S* 
along this path. The Hamiltonian associated with C is 



H=tf Q 



Pp) 



PQ K Q+PP K P 



(4) 
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TL is often called the "auxiliary" Hamiltonian, to avoid 
confusion with g. The points (Qfp,Pfp,Pq = 0, Pp = 
0) are fixed points of the auxiliary dynamical system 
formed from Tl. Let's note that there are infinitely 
many optimal trajectories from attractor to the sad- 
dle - they all differ by their initial momenta, but 
only paths that have zero energy give the stationary 
p(Q, P). One way to see this is to substitute p(Q, P) = 
p'(Q, P) exp {-S(Q, P)/D) into the stationary FPE, and 

obtain 22] (f§) 2 + (ff) 2 + ™K Q + $K P = 0. This 

is just a Hamilton- Jacobi equation, TL(Q, P, §p) = 0, 
and it describes dynamics on the "H = manifold. Now, 
since K(Qpp, Pfp) = 0, the momenta at all the FPs 
must also be zero. 

In summary, in the limit of weak noise, the rate of 
escape W ~ exp (— R/D) where R — S* - the action 
along the trajectory which connects the attractor to the 
saddle and satisfies H = 0. 



A. Analytical Limits 

Both limits discussed below are physically analogous 
to the two types of limits considered by Kramers [l[ : an 
overdamped limit and a low damping limit. One major 
difference is that our Hamiltonian g(Q,P) does not have 
an explicit separation of kinetic and potential energies. 



1. Overdamped regime in the vicinity of a bifurcation. 

At either \J~i^g or J (3g, the saddle and one of the 

attracting nodes of K merge together at (Qb,Pb) in a 
saddle-node bifurcation. As \f(5 approaches s/Pb, a sad- 
dle and a node approach each other, the attractive eigen- 
vector of the saddle, and one of the attracting eigenvec- 
tors of the node align parallel to each other, and the cor- 
responding attracting eigenvalues of both FP approach 
the same negative value. Also, the repulsive eigenvec- 
tor of the saddle and the other attracting eigenvector of 
the node become equal and the corresponding eigenval- 
ues both tend to zero. This process is depicted in the 
cartoon in Fig. [3J 

-+ 

Cartoon of K-dynamics near a bifurcation 




■ saddle 

• (attracting) node 



FIG. 3: Cartoon of the A-flow near a saddle-node bifurcation. 



Thus, close to a bifurcation, there develops a soft mode - 
a narrow region connecting the saddle with an attractor 

along which the A^-motion is slow. Due to this slow- 
ing down, the soft mode forms a path of least resistance 
along which a large noise-induced fluctuation away from 
the attractor is most likely to take place. To analyze 
this system we rewrite K in terms of parameter fj and 
variables (/ = Q — Qb,s = P — Pb)] it turns out to 
be unnecessary to make a linear transformation to eigen- 
coordinates of K at (Qb, Pb , VPs) ■ In terms of these 
new coordinates, K transforms to the following form: 

/ = -2fi- 1 (/ -a B s) + P B f 



- 2Q B sf + 3P B s 2 + fs + s 3 + n Q 



= ±fj-3Q B f -2P B fs 



) B s 2 - f - s 2 f + n P 



Here the + sign applies for the low field (high amplitude) 
bifurcation and 



a B = ± 2>/l - 3f2- 2 )/3. 



(7) 



Notice that the position of fixed points sets characteristic 
scaling of both s and / variables to be oc ry 1 / 2 . From this 
we see that s is a slow variable while / is fast - on its time 
scale, s appears approximately frozen. To lowest order 
in i), we treat s as completely frozen, in which case / 
relaxes to obs + corrections of order 77. In such adiabatic 
limit, the dynamics of s is then given by 



where 



±77- b(tt)s 2 + n P 



b(Q) = '3Q B a 2 B + 2P B a B + Qb 



(8) 



(9) 



(b > for low field bifurcation and b < for high field 
bifurcation). The potential for this slow coordinate is 



U = TV S + 77 s3 



(10) 



w hile th e distance between the saddle and the attractor is 
y/fj/\b\, thus the potential barrier for the slow coordinate 
is 



AU 



3\b\ x / 2 (ny 



(11) 



This is the well known £ = 3/2 scaling law. A 1- 
dimcnsional FPE for s can be written, a stationary 
solution for p(s, f), fl) obtained, and the escape rate 
calculated. This quantity will have the form W = 
W'(fi) exp (— AU/D). Such adiabatic approach is valid 
only sufficiently close to the bifurcation point. When 
fl <C 1, "sufficiently close" means that fj <C f2 -3 . In this 
regime, 



R : 



^fi" 1 / 2 ?? 3 / 2 low field bifurcation 



f2 1 r) 3 / 2 high field bifurcation 



(12) 



A more systematic treatment of the close-to-a-bifurcation 
regime which will allow us to extract corrections to the 
adiabatic limit and derive the fj f2 -3 criterion will be 
considered in Section III. A 



(5) 
(6) 
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2. Low dissipation regime. 

When effective friction, f2 _1 is low, the dynamics of 
K is approximately Hamiltonian: the time scale for the 
decay of energy is much larger then the time scale to 
make one cycle on the contour of approximately fixed 
energy. This separation of time scales breaks down very 
close to a bifurcation point, so as long as SI -1 ^ 0, 
there will always be a small fraction of the hysteresis 
displaying a £ = 3/2 scaling. Aside from this very nar- 
row region, the separation of time scales allows one to 
turn a 2-variable FPE into a 1-variable FPE for the 
diffusion of energy, by averaging out the fast part of 
the dynamics. The noise strength D is proportional to 
f2 -2 , so the FPE in the stationary regime can be rewrit- 
ten as {g,p} = -Q _1 V ■ (fp + D'D,~ l Vpj. Note that 

when f2 _1 = 0, any function of the form p{g{Q, P)) will 
be a solution to this FPE. Therefore, we can expand 
p(Q, P) = (g(Q, P)) + fry 1 ) (Q, P) + 0(n- 2 ) and 
integrate over a contour of constant g(Q, P). We will re- 
fer to such theory as " quasi-Hamiltonian" theory. After 
some algebra, this leads to: 



(13) 



where B{E) = f f g , Q P)=E dQ dP and D(E) = 
J J g ^Q P j =E \7 2 g(Q,P)dQdP. An alternative approach 
followed by [13] is to calculate an average rate of energy 
decay due to friction and an effective diffusion coefficient 
for the energy drift due to noise. These would turn out 
to be precisely the B(E) and D(E) respectively, hence 
Eq. (fl3|) is the FPE for the energy. Returning back from 
D' to D, the activation rate is W oc exp (-R/D) where 



R 



-fr 



B(E') 
D(E') 



dE' 



(14) 



Here E a is the energy of the attractor, and E s is the 
energy of the saddle. Close to bifurcation points, but not 
close enough for the overdamped theory to be applicable, 



R : 



-1^3/2 



low field bifurcation [18] 
high field bifurcation [24J 



(15) 



Comparisson of Eqn. (fl"5|) with Eqn. p2[) reveals that 
there must be a crossover from a £ = 1 regime to £ = 3/2 

regime as a/p approaches 



B. Power law scalings of escape barriers 

When SI -1 ^ 0, sufficiently close to the saddle- node 
bifurcation, the dynamics near the attractor and the sad- 
dle is generically over-damped [2!| . Hence, the soft-mode 
geometric picture offered in Section II. A. 1 is generic and 
£ = 3/2 is a genetic property in a system with finite 



damping and close enough to the saddle-node bifurca- 
tion. Further from the bifurcation, the over-damped soft 
mode does not exist. We have seen that for sufficiently 
low f2 _1 , the overdamped £ = 3/2 crosses over to £ = 1 



(close to 



or remains £ = 3/2 (close to 



As 



recently proved by Dykman, Schwartz and Shapiro [271 ] 
(DSS), these are in fact also generic scaling laws in the 
under-damped vicinity of a saddle-node bifurcation. DSS 
considered two scenarios which may happen when f2 _1 
equals zero. As f) — > 0, the center and the saddle may 
merge, or they may remain separate. DSS named the 
first type of a Hamiltonian bifurcation "local" and the 
latter " non- local" . The two situations dictate the form 
of the Hamiltonian g. 

Let's switch to coordinates (q,p) such that when 
SI -1 = 0, both the center and the saddle fixed points 
lie on the p — axis. In the case of a local Hamiltonian 
bifurcation, in the vicinity of (q — 0,p = 0, fj = 0), g 
is given by gL(q,P,v) = \p 2 + U cub {q) where U cub (q) = 
— ig 3 + fjq. The center lies inside a small homoclinic 
loop at the energy of the saddle. DSS proved that given 
this form of the Hamiltonian g, addition of noise and 
damping leads to R{f)) = A(Q)rj 3 / 2 (in the under-damped 
regime). But we know that as as the system approaches 
the bifurcation even closer, it becomes over-damped, a 
soft-mode sets up and the method of analysis presented 
in Section II. A. 1 leads to R(fj) = _B(fi)ry 3/2 . Therefore, 
assuming smooth dependence upon fj, it must be that 
A(Cl) = B(tt). In other words, the £ = 3/2 power law 
in this situation holds beyond where it was expected to 
hold due to a soft mode theory! In the Duffing system, 

a nonlocal bifurcation takes place at J f3 B . In the nu- 
merical experiments described below, we found that the 
£ = 3/2 power law holds near y /3 B over a finite range of 

fj even as S! -1 — > 0. This £ = 3/2 law at zero fi -1 has 
been predicted by 24], but the new work of DSS gives 
it a new explanation. We elaborate upon this in Section 
III.B. 

In the case of a nonlocal Hamiltonian bifurcation, g 
is given by g NL (q,p,fj) = \p 2 + fjU(q), where U(q) is 
independent of fj and has a local minimum and maxi- 
mum. The homoclinic loop of such systems has the prop- 
erty that only its size in the p-direction shrinks to zero 
as i) — * 0, while the distance between the center and 
the saddle inside the loop remains finite. Introduction 
of damping into these types of systems will not change 
the separation between the saddle and the attractor un- 
til the system has come close enough to the bifurcation: 
fj ~ fl~ 2 . DSS proved that given this form of the Hamil- 
tonian g, addition of noise and damping leads R{fj) cx fj 
(in the underdamped regime). Therefore, at finite f2 _1 
we expect a £ = 1 scaling moderately close to the bifur- 
cation, but further away then the region of over-damped 
£ = 3/2 scaling. In the Duffing system, a nonlocal bifur- 
cation takes place at \ (3 l B . It is important to stress that 
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all saddle-node bifurcations at finite f2 _1 are local in the 
sense that the attractor and the saddle merge as f\ — > 0. 
The terms "local" or "nonlocal" only make sense at zero 
f2 _1 . We must also mention that a nonlocal bifurcations 
may take place in many forms. For example, in a Duff- 
ing oscillator, the homoclinic loop has a horse-shoc shape 
at Q^ 1 = which becomes thinner while remaining ex- 
tended as fj — > 0. DSS found a change of variables which 
maps from g to g^L- 



III. NUMERICAL METHOD AND RESULTS 



regime, finding a trajectory between an attractor and the 
saddle is difficult only very close to bifurcation points, 
but in the large- regime, finding such a trajectory is 
difficult over a larger part of the hysteresis. 

We will now summarize results of the computation of 
S* , scaling laws, and crossovers. It is sometimes easier to 
present results in terms of x = — ==^ — p== - "reduced" fj. 

For each fl~ 1 and each type of transition, we scanned the 
hysteretic region at multiple values of x and for each x we 
followed the "shooting" procedure, and where necessary, 
the bracketing method. 



Numerical calculations of escape rate are based on di- 
rect evaluation of S* . The algorithm for doing so is as fol- 
lows. (1) Choose parameters and ^//3. (2) Calculate 
repulsive eigenvectors (i.e. those whose eigenvalues have 
positive real part) oftheFP (Qattr, P a ttr,PQ — 0,pp = 0) 
of the auxiliary system. These eigenvectors span the tan- 
gent plane to the unstable manifold of this attractor. (3) 
Position a locus of initial conditions on a small circle 
that lies on this tangent plane and centered around this 
FP and evolve each of these initial conditions according 
Hamiltonian equations with Hamiltonian Ti from Eqn. 
0. The angle around this circle, <p, serves as a param- 
eter that enumerates a trajectory. (4) From this set of 
trajectories, find the one which leads into the saddle, with 
a special trajectory parameter ip* . Because the method 
involves "shooting" many trajectories to see which one 
hits the saddle, this method is occasionally called "the 
shooting method". (5) Evaluate the action along this 
trajectory. For each fi , this procedure would be re- 
peated for many values of y/]3 between and *J~0%- 

In practice, step (4) may require several rounds of 
bracketing. In each round, for each trajectory, we cal- 
culate the distance between the saddle and the point of 
intersection of that trajectory with the separatrix, mea- 
sured along the separatrix. We identify a p for which 
the trajectory came closest to the saddle and then evolve 
another set of trajectories in a small range of trajectory 
parameters around this p. After more and more rounds 
of such procedure, we would obtain a trajectory that hits 
the separatrix closer and closer to the saddle. In general, 
finding a trajectory that connects an attractor with the 
saddle is difficult because it lies on the intersection of the 
unstable manifold of the attractor and the stable mani- 
fold of the saddle in the 4D auxiliary space [13], [13] ■ A 
small deviation of p on opposite side of (p* will produce 
trajectories that lie on opposite sides of this intersection 
and will diverge from each other. This effect is most pro- 
nounced closer to bifurcations because the motion along 
the optimal trajectory becomes much slower then the di- 
verging motion away from it - such is the structure of 
eigenvalues of both attractor and saddle close to the bi- 
furcation. This polarity takes place only close to bifurca- 
tion points, but it persists further away from bifurcation 
points when f2 _1 is larger. Therefore, in the small-f^ 1 



A. Up — > Down transitions 

1. Numerical results 

A plot of S* /Vl -1 versus x, Figs. [4] [8] reveals several 
features. At low SI -1 two regimes clearly stand out: one 
with £ 3/2, and one with £ rs 1. Further from the bi- 
furcation, the £ = 1 scaling law breaks down and S'*/0 _1 
versus x displays a non-power-law behavior. For various 
values of fi -1 , the data in the "3/2 regime" fits well to 

a) S>„' 
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FIG. 4: Main plots: S^/fT 1 (dots) and fits (lines) for up 
— > down transitions. The fit to the £ = 3/2 regime was made 
by analyzing In (S 1 */^ 1 ) vs. In x and a fit to the £ = 1 regime 
was made by analyzing S* /QT 1 vs. x. Inserts: S* (i)/Sl _1 
(dots) and theory (lines) based on Eqn. [12] for the £ = 3/2 
regime and Eqn. [13] for the £ = 1 regime. 
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a power law with £ which is up to 3% below the theo- 
retical value of 3/2. In general, precise extraction of the 
exponent 3/2 requires approaching very close to the bi- 
furcation. It appears that the discrepancy was highest 
in those cases when the bifurcation was approached less 
closely (relative to the crossover point) or less reliably. 
The reliability of the approach is hampered at larger fi -1 
by the strong divergence of trajectories explained in the 
previous section. A crossover between the £ = 3/2 and 
£ = 1 regime is defined to be such x at which fits to the 
respective regions intersect. Collecting data for S* vs. 
x for various we were able to map out the entire 
bi-stability region for locations of various scaling regimes 
and crossovers between them. The lowest at which 
computations were made was 10 -3 . The product of this 
work is depicted in Fig. [5) It is important to stress that 
in Fig. [5] for up — > down transitions and in Fig. [9] for 
down — > up transitions, the mapping was performed by 
scanning y/]3 while holding fixed. 




FIG. 5: Regions of different scaling behaviors for up — > down 
transitions obtained by scanning y 7 /? at fixed at multiple 
values of Region I - non-power-law regime. Region II - 

£ = 1 regime and region III - £ = 3/2 regime. The crossover 
between the £ = 1 and a non-power law regime was defined 
as the point at which the fit to £ = 1 regime came closest to 
the In (Q'/Q.- 1 ) vs. In x - see Fig. [4] 

Note that as f2 _1 — > 0, the linear scaling is found between 
x = and some finite x. At small, but non-zero 
the region of £ = 3/2 subsequently grows. These results 
are in accord with the theorem of DSS which states that 
close to a bifurcation which at f2 _1 = becomes nonlo- 
cal, in the underdamped regime at non-zero the £ 
is expected to be 1. Since the underdamped portion of 
the hysteresis extends all the way to the bifurcation as 
— ► 0, we expect the linear scal ing to extend all the 
way to the bifurcation in this limit [331 ] . 

Clearly, the boundaries between regions I, II and III 
are purely conventional - they have been defined in some 
convenient way, but do not correspond to anything specif- 
ically physical. To emphasize this point, it helps to con- 
struct a 3-dimensional plot of S* /Q^ 1 versus x and SI -1 



- Fig. [6j This surface is smooth and large parts of it 



0. 




FIG. 6: S'/Q- 1 versus x and Q.- 1 up to fT 1 = 0.175 which 
is about 30% of the hysteresis. Fig. [4] are two slices of this 
surface. 

happen to fit well to power laws. The regions when this 
is possible appear to be much bigger then the regions at 
which the £ = 1 and £ = 3/2 scaling was supposed to 
be applicable on the basis of approximations discussed in 
Section II. A natural question is why this is so! To shed 
more light on this issue, let's digress on the crossover 
from the £ = 3/2 to £ = 1 in more detail. 

2. Discussion of the £ = 3/2 to £ = 1 crossover 

Notice that the crossover between the £ = 3/2 and 
£ = 1 regime, as it appears in Fig. [4] was defined as the 
point of intersection of the fits to the respective regimes. 
Thus, to predict this crossover we simply equate the 
asymptotic in the £ = 3/2 regime, R = ^-VL~ 1 l 2 ff/ 2 
to the asymptotic in the £ = 1 regime R — 2fl~ 1 fj 
(this is most accurate as SI -1 — > 0). The result is 

fj = which corresponds to \f]3 — \j^ B + 

(^is itself a function of fl"*). The cross-over pre- 

dieted this way is depicted by dashed line in Fig. 5. The 
fact that this method of predicting a crossover works well 
tells us that the asymptotic expressions are correct. How- 
ever, we could imagine that the width of the crossover 
region was very large, in which case the crossover de- 
fined by equating the fits and predicted by equating the 
asymptotics is not physically meaningful. The width of 
the crossover region is a more physically interesting quan- 
tity. Alternatively we may also define the following two 
crossovers: a "lower crossover" (LC) defined as a charac- 
teristic value of \/]3 (or fj or x) above which the £ = 3/2 
regime begins to fail and an "upper crossover" (UC) de- 
fined as a characteristic value of \fp (or fj or x) below 
which the £ = 1 regime begins to fail. The width of 
the crossover region is the difference between these two 
values. 
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Physically, we expect the LC to take place where the 
soft-mode picture fails. Therefore, a good characteristic 
estimate of this will be such fj at which the real and 
imaginary parts of eigenvalues become equal. This can 
be shown to happen for fj = £l~ 3 /2 - the scaling is cubic. 



To predict the LC more accurately let's turn back to the 
type of analysis described in section II. A. 1. Let's re- 
write Eqns. (5)-(6) for the non-adiabatic part of the fast 
variable, A = /-ags. Near low- field bifurcation, these 
will have the form 



A = -a B fj-2Q l A + ais 2 + a 2 As + a 3 A 2 + a 4 s 3 + a 5 s 2 A + a 6 sA 2 + a 7 A 3 + n x - a B n y (16) 
s = fj + hs 2 + b 2 As + b 3 A 2 + b 4 s 3 + b 5 As 2 +b 6 A 2 s + b 7 A 3 + n y (17) 

I 



All coefficients a; and bt are functions of Or 1 . The posi- 
tion of fixed points sets the characteristic scaling of A to 
be fj and of s to be r) 1 / 2 (generic feature of saddle-node 
bifurcations), so we define new variables: s — fj l l 2 s' and 
A = fjA' . Moreover, the ^-dependencies of the coeffi- 
cients ai and bi near low-field bifurcation are such that 

I 

d<p 

lw 

do 



where r = Sl~ x T . With the help of the definition of a 
and Eq. ([7]) we note that the zeroth-order part of the 
fast variable o^s is given by f^^V^^a and the cor- 
rection due to non-adiabaticity of the fast variable is 
A = ijn 2 (j). We also see from Eqs. (ll~8l) - (fT9|) that when 
(j, -C 1, both a and <fi are 0(1). Hence comparison of 
the two terms says that when /i <C 1, non-adiabatic 
part/adiabatic part ~ {f)Vl 2 ) / (f] 1 / 2 ^}/ 2 ) = [i x l 2 -> 
as [i — ► 0. In other words, in this regime the non- 
adiabatic part becomes unimportant and the problem 
reduces to 1 dimension defined by the noisy dynam- 
ics of the slow variable. To do this reduction explic- 
itly, we first notice from Eqs. (fT8"|) - (fll?)) that at lowest 
order, the variable a can be considered constant while 
(f>(cr) ~ \o 2 — \ + non-adiabatic corrections. This 4>{<j) 
is then substituted into Eq. (fl~9)) . The resulting expres- 
sion can only be considered up to 0(/j,) since the higher- 
order terms must also include the non-adiabatic correc- 
tions to 4> to be complete [note that if the /i 1 / 2 -terms in 
Eq. (|19p contained any <p terms, the complete expression 
within the adiabatic approximation would only extend 
to C^/i 1 / 2 ), not O(n)]. The barrier height in the resul- 
tant 1-dimensional potential is All — ^^[i 1 ! 2 + \f2\j?l 2 . 
Again, the correction is within the adiabatic approxima- 
tion; a non-adiabatic correction would be higher order 
in fi. Taking into account that the N„ has a diffusion 



when fl^ 1 is small compared to 1, further re-scaling of 
s' and A' as s' = Qr x l 2 o and A' = fl 2 (f> yields a single 
dimensionless parameter \i = fjQ 3 ; such reduction is ap- 
proximate, but becomes more accurate for smaller f^ 1 
[34j |. In terms of n, a and <j> we have 



(18) 
(19) 



constant D a = ^-D, we have 

r = ^n-i/ 2 7 f/ 2 + V2ft 5 /V /2 (20) 

o 

[compare with Eq. (|12j) ] . We see again that the correction 
term grows to be dominant when /i ~ 1 (i.e. when fj ~ 
Q,~ 3 ) , but this is precisely when the non-adiabatic part 
becomes comparable to the adiabatic part, and it turns 
out to be meaningless to treat the problem as effectively 
1-dimensional. Our analysis taught us that in the regime 
when fj <C f^~ 3 , the activation barrier R will scale as fj 3 ^ 2 
to a good approximation. Outside the regime of 77 <C £l~ 3 
the current analysis only allows us to conclude that if £ 
remains to be 3/2, the physics behind this scaling is not 
a physics of a 1-dimensional over-damped soft mode; the 
escape problem is not 1-dimensional for fj 3> £l~ 3 . 

To answer the question of what happens outside of the 
ry 3> Q~ 3 regime we chose to compute the most likely 
escape path numerically and calculate S* on that path. 
We looked at the difference between the zeroth-order the- 
ory, i?o = -M^f^ 1 / 2 ?^ 3 / 2 and the numerically-computed 
S* . This was repeated for several values of f2 and in each 
case, a fit to the function oc ff was made. Although this 
represents computer calculations, the data is somewhat 
noisy due to the difficulty of hitting the saddle. Values 
of p between 2.15 and 2.86 were found. The predicted 
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value of p = 5/2 in the correction factor lies in the mid- 
dle of this range. However, the crossover LC, defined as 
an intersection of some small fraction r of Rq with a fit 
to i?i scales linearly with O -1 , not cubically. The fj^c so 
obtained can be fit to a power law with exponent 1.2 for 
r = 0.1 and and 1.1 for r = 0.2 - both of these exponents 
are far from 3. 



The question of why the £ = 3/2 scaling holds over 
a much larger region of parameter space then expected 
based on the 1-dimensional soft-mode picture remains 
unanswered, and forms a good challenge problem for a 
future work. One potential hypothesis is that in the re- 
gion when the 1-dimensional soft-mode picture does not 
hold true, the action along the MPEP actually grows 
only over a small portion of the MPEP where the mo- 
tion is essentially 1-dimensional. We provide plot of the 
action versus the distance along the trajectory, S(£) and 
the derivative dSjdt for a particular value of parameters 
deep in region III - Fig. [7j The derivative plot clearly 
shows that action grows fastest in the last winding of the 
trajectory, but this is not sufficient to conclude that the 
dynamics in that part of the trajectory can somehow be 
described as effectively 1-dimensional. 



Similarly to the LC, the UC can be defined as a char- 
acteristic value of the driving field below which the lin- 
ear scaling begins to break down. To predict this point, 
we must again compare the lowest and one higher order 
terms in the theory for R(fj) that applies in the £ = 1 
regime. The effect of a non-zero fi -1 is clearly seen 
in Fig. [5] The lowest-order in dissipation theory (quasi- 
Hamiltonian) indeed predicts linear scaling for low values 
of rj (or x), which according to Eq. (fTS"|) is given by 2£l~ 1 fj. 
The region of linear scaling persists for non-zero fi , but 
due to the existence of the 3/2-scaling region, this linear 
part has been "pushed over" to larger rj and to lower 
values then those predicted by the zeroth-order theory. 
A work by Chinarov' et. al. [25j worked out dissipative 
corrections to the quasi-Hamiltonian theory to next order 
in f2 _1 . Translating their results to our notation reads: 
R = tt^ 1 (2fj - 7rSl _1 ) (compare with Eq.[T5]). The factor 
7ril~ 2 is just that correction which takes this lowering into 
account. We expect the linear scaling to start breaking 
down when the correction reaches some fraction r of the 
zeroth-order term, i.e. when fjuc — ~ 1 f° r 

low fi . To extract fjuc from the data we subtracted the 
zeroth-order part from S* /tt -1 . The resulting function 
grows for low 77, slows down for intermediate fj (but not 
completely because the exact slope in the linear regime 
is slightly different from the zeroth-order result) and in- 
creases for yet larger fj past the linear regime (see Fig. 6). 
The average value in the intermediate regime is taken to 
be the offset factor. Equating this correction to the ze- 
roth order term yields fjuc which does depend linearly 

onrr 1 . 
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FIG. 7: S(l) and its derivative for particular parameters in 
region III as well as the plot of the MPEP (dashed curve 
denotes separatrix). Here = 0.15 and \fP = 10% of 

the hysteresis. These parameters correspond to /j, w 5; the 
soft-mode picture is expected to hold only for 



B. Down — > Up transitions 

Analogously to the up — > down transitions, the en- 
tire bi-stability region was mapped out. The product 
of this mapping is shown in Fig. [5] The lowest fi" 1 at 
which computations were made was 0.005. Fuzzy re- 
gions correspond to those sections of the - O^ 1 plane 
which proved exceptionally difficult to map reliably due 
to high divergence of trajectories. A portion of the sur- 
face S* /VL^ 1 {x, is presented in Fig.fTOl 

The low amplitude basin does not exhibit nonlocal ge- 
ometry upon approach to the bifurcation at zero f^ 1 , 
and correspondingly, the linear scaling is not found in 
down — > up transitions at any f2 _1 (see Fig. Hip , which 
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S>n' 




FIG. 8: R* '(x) /f2 _1 on a linear scale. The thick black line 
represents the zeroth-order theory, Eq. (I14|l . The other curves 
represent numerical calculations of S* - squares: SI -1 = 0.03, 
circles: fi _1 = 0.06. 



s /n- 




0.0 



FIG. 10: S'/Q" 1 versus x and fT 1 up to fT 1 = 0.22 which 
is about 33% of the hysteresis. Fig. [4] are two slices of this 
surface. 




FIG. 9: Regions of different scaling behaviors for down — > up 
transitions obtained by scanning y 7 /? at fixed Or 1 at multi- 
ple values of f2 _1 . Region I corresponds to a non-power-law 
regime. Region II corresponds to the £ « 1.3 regime. Region 
III corresponds to the £ = 3/2 regime. 

is in accord with the DSS theorem. The £ = 3/2 scaling 
gives way to a different scaling relationship further from 
the bifurcation, with an exponent £ ~ 1.3. We interpret 
this £ 1.3 exponent as a power-law interpolation. 

Note that region III in Fig. [5] (the £ = 3/2 region) 
extends all the way to zero f2 _1 - this is unlike section 
II, where the £ = 3/2 region disappears as SI -1 — > 0. 
These results are in accord with the theorem of DSS 
which states that close to a bifurcation which is local 
at il^ 1 = 0, in the underdamped regime at non-zero f2 _1 
the £ is expected to be 3/2. Since the underdamped 
portion of the hysteresis extends all the way to the bifur- 
cation as £1 _1 — *• 0, we expect the 3/2 scaling to extend 
all the way to the bifurcation in this limit. Had we as- 
sumed that £ = 3/2 due to the over-damped soft mode, 
we would follow the nearly-adiabatic analysis of the pre- 
vious section and we would find that for low fl _1 there 
exists another dimensionless parameter v = fj 1 / 2 ^ 2 (35|. 



When i/<l such analysis predicts that £ = 3/2. When 
v ~ 1, the reduction to 1-d due to adiabaticity becomes 
impossible and prediction based on this method can not 
be made. The fact that £ remains 3/2 beyond v <C 1, 
as we see from numerics, indicates that the 3/2 scaling 
is not a consequence of reduction to an over-damped 1-d 
soft-mode. 

In the previous section we offered three definitions of 
the 3/2 — > 1 crossover. In the case of up — > down tran- 
sitions, we don't have theory for the £ ss 1.3 scaling 
regime, so the crossover from the 3/2 scaling can be de- 
fined here only as a point at which the 3/2 scaling be- 
gins to break. Since the 3/2 scaling is predicted by the 
quasi-Hamiltonian theory, we look at next order correc- 
tions due to Chinarov et. al 25]. Translating the results 
of that work to our notation, close to the bifurcation, 
R w gTTjft-V/ 2 + 1.06g4?n- 3 llnO" 1 ! rfl 2 (compare 
with Eq. 1 1 5|) . Evidently, the correction is to simply shift 
R(fj) by a constant offset (see Fig. [TT] (b) where this ef- 
fect is clearly pronounced throughout the entire R(fj), not 
just close to the bifurcation). Because the correction is 
so small for small f2 _1 , its effect doesn't become notice- 
able until a rather large f2 _1 - equating the zeroth-order 
term to the correction term predicts this value of fl^ 1 to 
be w 1.5. Qualitatively, this delayed effect is observed in 
numerics (see Fig. [SJ, but the crossover doesn't seem to 
depend on S! -1 significantly until fi" 1 ss fi" 1 /2 w 0.28. 



IV. CONCLUSION 

We considered the escape rate W cx exp (-R/D) from 
one basin of a non-linear oscillator into another. We cal- 
culated the leading term in the effective barrier R in the 
limit of small noise: In W = -i?(rr\ yf/3)/D + O(D ). 
For any given ft -1 , this R may exhibit several scaling 
behaviors versus | — \/ Pb\ where V/^b is a value of 
the driving field at which a saddle-node bifurcation takes 
place. We used numerical shooting method to calculate 
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a ) S/h< 




0.001 0.01 0.1 1 



b) S> n - 




FIG. 11: Main plots: S*(x)/Q~ 1 (dots) and fits (lines) for 
down — * up transitions. Fits to the £ = 3/2 and £ ~ 1.3 
regimes was made by analyzing In (S* /Q." 1 ) vs. lnx. The 
region where the simple shooting method was unreliable for 
hitting the saddle (due to high divergence of trajectories), a 
bracketing method was used, the results of which are denoted 
by squares. Inserts: S* (x)/Q~ 1 (dots) and theory (lines) 
based on Eqn. [14] for the £ = 1 regime; note that this quasi- 
Hamiltonian theory predicts the £ — 3/2 scaling, as was first 
pointed out by [24[ and in light of the recent work by [23] this 
fact has been explained to be a result of the fact that this 
bifurcation is local at = 0. 

most probable escape paths (MPEPs) from one basin to 
another. R is the action along this MPEP between an 
attractor inside one of the basins and the saddle that 
lies on the separatrix between the two basins. This tech- 



nique allowed us to map out the v 7 /? - parameter 
space for locations of various scaling behaviors of R ver- 
sus — sffis • The result of this mapping is shown in 
Fig. for up — ► down transitions and in Fig. [5] for down 
— > up transitions. 

We confirmed predictions of the recent work by Dyk- 
man, Schwartz and Shairo (DSS) [13, which classifies 
scaling laws based on the structure of the homoclinic tra- 
jectory at the energy of the saddle zero = - "lo- 
cal" or " nonlocal" . We found that for those bifurcations 
which become nonlocal when SI -1 = 0, at non-zero O x 
there exists a regime with R cx | ^f]3 — \[]5b | > as predicted 
by DSS. Such situation occurs for up — * down transi- 
tions close to the low-field bifurcation. Even closer to 
the low-field amplitude bifurcation we found a crossover 

into the R cx — VPs | 3 ^ 2 regime. We have defined 
this crossover in three different ways. Following one of 
the definitions we discovered that the 3/2-scaling breaks 
down much further from the bifurcation then where it 
would be expected to break down if the physics of this 
scaling was due to an over-damped soft mode. 

We also found that for those bifurcations which are lo- 
cal when Q^ 1 = 0, at non-zero il^ 1 there exists a regime 

with R oc y/]3 — V '/3b , which holds over a finite po- 
tion of a hysteresis even as — > 0. Such situations take 
place for down — > up transitions close to the high-field 
bifurcation (this has been predicted earlier [Hj]). This 
3/2 scaling again takes place much further from the bi- 
furcation then were it would be expected to hold as a 
result of over-damped soft mode, but for the down — > 
up transitions this fact is somewhat less novel because 
it has been explicitly predicted to hold in the limit of 
— ► a long time ago by [24[ and also by recent work 
of DSS as explained in Section II. B. Further from the 
high-field bifurcation, we found another crossover into 
a regime 

R oc vWif where £ w 1.3. A crossover 

into power law regime was not expected on any analytical 
grounds, so we interpret this power law as an interpola- 
tion. 
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